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ABSTRACT 

We present ground-based broad-band photometry of two transits in the WASP-44 
planetary system obtained simuhaneously through four optical (Sloan g' , r' , i' , z') 
and three near-infrared (NIR; J, H, K) filters. We achieved low scatters of 1-2 mmag 
per observation in the optical bands with a cadence of « 48 s, but the NIR-band light 
curves present much greater scatter. We also observed another transit of WASP-44 
b by using a Gunn-r filter and telescope defocussing, with a scatter of 0.37 mmag 
per point and an observing cadence around 135 s. We used these data to improve 
measurements of the time of mid-transit and the physical properties of the system. In 
particular, we improved the radius measurements of the star and planet by factors of 3 
and 4, respectively. We find that the radius of WASP-44 b is 1.002 ±0.033 ±0.018 i?jup 
(statistical and systematic errors, respectively), which is slightly smaller than previ- 
ously thought and differs from that expected for a core- free planet. In addition, with 
the help of a synthetic spectrum, we investigated the theoretically-predicted variation 
of the planetary radius as a function of wavelength, covering the range 370-2440 nm. 
We can rule out extreme variations at optical wavelengths, but unfortunately our 
data are not precise enough (especially in the NIR bands) to differentiate between the 
theoretical spectrum and a radius which does not change with wavelength. 
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1 INTRODUCTION 

The transit method is not only an excellent technique to 
detect extrasolar planets, but also provide unrivalled access 
to their physical parameters. The geometry of a transiting 
extrasolar planet (TEP) system permits measurement of 
the mass, radius and density of a planet. These quantities 
require some input from stellar evolutionary theory, but the 
planetary surface gravity can be measured directly. A TEP 
system is therefore a potential ho rn of plenty of physical 
information on planet a ry systems ( Seager &: Sasselovl 2000l: 



2003; Sudarsky et al.' 
Holman & Murray 2005; 



Sudarsky et all l2000l: iBrown et al.| 1200 ll : iHubbard et al 
200l| : ISeagcr fc Mallen-Ornelas' 

2OO3I: [Char bonncau ot al-, .2005; 

Winn et al.ir2005l : ISouthworth et al]i2007l ). 

Currently, most of t he TEPs discover e d by ground- 
based transit surveys (e.g. iHellier et al.]|2012l : ISmallev et al.l 
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1201^; 'Bakos et al.' '2012^; 'Hartman et al.' '2012"; 'Pene v et all 
2013; Bryan et al. 2012; Bcatty ct al. 2012; Sivcr d et all 
[20121 ) are close-in hot Jupiters, because the detection prob- 
ability for such planets is much greater than for those 
which are smaller or on wider orbits. Hot Jupiters are also 
relatively straightforward to detect via the radial velocity 
method, as they induce comparatively large velocity vari- 
ations in their h ost stars. The existence of hot Jupiters 
was unexpected l|Mayor and Queloj|l995l ). and led to the 
birth of new formation and mi grati on theories (see the re - 
views of iD'Angelo e FaT '2010' and iLubow and Idal I2OI0I '). 
The bias of the ground-based surveys in favour of this class 
of planets is clear, whereas the Kepler satellite has already 
found hundreds of smaller and long er-period planet cand i- 
dates from only 16 months of data jBorucki et afllio 11 al ibi: 
iBatalha et~alll2012l '). 

Another fascinating opportunity offered by TEPs is the 
possibility to probe their atmospheres. Due to atomic and 
molecular absorption, the analysis of the spectrum of the 
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parent stars during transit (primary eclipse) is an excellent 
way to probe their atmospheric composition. The existence 
of two different classes of hot Jupiters (pM and pL) has been 
suggested according to the incident stellar flux received from 
the parent star and to the expected amount of absorbing 
substances, such as gaseous titanium o xide (TiO) and vana - 
dium oxide (VO), in their atmospheres (jFortnev et al.|[2008h . 
According to theoretical models, these two oxidized elements 
should be extremely strong absorbers for the pM class be- 
tween 450 nm and 700 nm, causing such planets to have a hot 
(~2000 K) stratosphere. An observable variation of radius 
with wavelength, due to the opacities of the se elements, is 
conse quently expected in the optical bands (|Burrows et al.l 
l2007h . The radius of a pM planet should be 3% lower at 
350-400 nm versus 500-700 nm. The radius variations for the 
colder pL-class planets (incident flux <10^ ergs~^ cm~^) are 
smaller, but a significant contribution to the opacity in the 
optical ba nd is expected from Na I at ~590 nm, and Kl a t 
-770 nm (|Fortnev et all |2008| . [2OI0I : iBurrows et al.ll201ol l. 
Transmission spectroscopy is essential for detecting atmo- 
spheric absorbers, and is complementary to observations at 
secondary eclipse which can help to determine the day-side 
atmospheric temperature structure. 

Accurate measurements of planet radii provide critical 
constraints for astrophysicists working on planet formation 
and evolution, because they provide an indication of the 
planet internal heat, structure and composition, as well as 
the heating mechanism for those which undergo strong tidal 
effects and stellar irradiati on, allowing the discriminatio n 
between competing theories jPollack et al.|[l99^ : lBosslll997^ . 
In particular, the most irradiated hot Jupiters are charac- 
terised by anomalously inflated radii for which t he expla- 
natio n remain s uncle a r jBodenheimer et al.ll2003l: Gu et al 



2004; Kiaudi '2005'; 'Fortnev et all 



200£ ; toobbs-Dixon & 



Demorv fc Seagejl201ll ) 



LiJ 

Ibgui fc Burrow j |2009| ; (Miller et al 



I - et ai.i 12OO7I; uacKson et ai. 
20081; lBa.tvgiii et al.l l2009l; 



Jackson et al 



20091; llbgui et al 



201 d; iLaughlin et all I2OIII ; iMiller fc FortnevI l20lil : 



Photometric and spectroscopic analyses of planet at- 
mospheres started with the use of optical, near-infrared 
(NIR) and infrared (IR) instruments aboard the Hub 



hie and Spitzer spa c e teles c opes fchar bonne au et al] 



Vidal-Madiar et all l2003l. l2004l: 



2002 



2ooe 



Tinetti et al.ll2007; 



2007 



Swain et al. 



Knutson ct al, 
2008; Pont ct al 



Richardso n et af 
2007a b; Richard son et a] 
2008; Beauheu et al 



20081 . 120101 . l201ll ; lGillon et al.ll201 j Tand led to the detection 



of Rayleigh scattering in a few hot Jupiter atmospheres, plus 
absorption lines associated with H2O, Na, CH4, TiO and VO 



20071; Sing et al.l 2008a bl; iLecavelier des Etangs et al..,200a 


Desert et al.1 20081; IZahnle et al.l 2009; Spiegel et al. 


200c 


Fossati et al.1 2010|; 'Tinetti et al.| 201GI; iDesert et al 


2011 


Wood et al.ll201ll; Gibson et al.ll201ll; ICrossfield et al. 


2012 


Crouzet et al.ll2012l). 



After a number of pioneering experiments with null 
results, ground-based spectroscopy also obtained some in- 
teresting results. Observations performed at the Hobby- 
Eberly and Subaru telescopes, in the spectral range 500- 
900 nm, detected Na absorption in the t ransmission spec- 
trum of HP 189733 b (jRedfield et al.ll2008l ) and HD 209458 b 
ijSnellen et al.ll2008l ). NIR spectroscopy of HD 209458 b at 
the VLT in the range 2291-2349 nm reported a significant 



wavelength shift in absorpt ion lines from carb on monoxide 
in the planet's atmosphere l|Snellen et al.ll201oh . 

High-resolution NIR spectroscopic measurements at the 
Keck telescopes over the wavelength range 2100-2400 nm, 
revealed an uppe r limit for CO a.bsorp tion in the atmosphere 
of HD 209458 b (jDeming et al.1 120051) . and that the super- 
Earth GJ 1214 b is a H-dominated planet (jCrossfield et al.l 
l201lh . GJ 1214 b was also studied through multi-object spec- 
troscopy from 0.61 to 0.85 fim, and in the J, H, and 
K atmospheric windows by using VLT/FORS and Magel- 
lan/MMIRS, the data being consiste nt with a featur eless 
transmission spectrum for the planet (|Bean et al.ll20lil ). 

A radius vari ation with wavelength was investigated 
for HD 209458 b bv lKnutson et al.1 ()2007al 'l using ten-colour 
HST photometry. Di fferent analyses of thes e data have given 
conflicting result s (Knutso n et~a\\ l2007al ; iBarmanI |2007| ; 
[Sing ct all l2008al ; [Southworthl I2OO8I ). An HST transmis- 
sion spectrum of HD 189733 b covering 270-570 nm showed 
a gradual increase of radius towards shorter wavelengths 
which was interpreted as a result of Rayleigh scatte ring from 
a high-altitude atmospheric haze (|Sing et al.ll2011al ). But ob- 
servations covering 550-1050 nm did not provide any indica - 
tion of the expected NaT or Kl features (jPont et al.ll2008l l. 
and those in the ranges 1082-1168 nm and 151 4-1693 nm did 
not de tect any variation in planetary radius (|Gibson et al.l 
l2012al ). 

Us ing the Gran Telescopio Canarias fGTC). fSing et aP 
(|2011bl ) claimed the flrst evidence for potassium in an extra- 
solar planet, from photometry of the XO-2 system in four 
narrow red-optical passbands, by detecting K l absorption 
at 766 .5 nm. This technique was also used by IColon et al.l 
(|2012i ) to study HD 80606 b, who found an unexplained large 
(« 4.2%) change in the apparent planetary radius between 
the wavelengths 769.9 and 777.4 nm. The differential spec- 
trophotometry technique was us ed to obt ain two transit light 
curves of XO-2b with GTC bv ISing et al.l (|2012l ). who de- 
tected signiflcant absorption in the planetary atmosphere 
in a 50- A bandpass centred on the Nal doublet. Instead, 
the presence of an Na -rich atm osphere in WASP-29b was 
recently ruled out bv iGibson et al.i (|2()12b), using Gemini- 
Sou th GMOS transit spectrop hotometry. 

ISouthworth et al.1 (|2012bl ) recently presented a study 
of possible radius variations of the hot Jupiter HAT-P-5 b, 
based on photometry obtained simultaneous in the u, g, r 
and / passbands. The authors detected a gradual increase of 
the radius between 450 nm and 850 nm, plus a substantially 
larger planetary radius at 350 nm. The latter phenomenon 
can be explained by systematic errors in the u-band pho- 
tometry, but is also consistent with Rayleigh scattering as 
in the case of H D 189733 b. A similar multi-band investiga- 
tion of HAT-P-8 jMancini et al.ll2013h suggests the presence 
of strong optical absorbers near the terminator of this TEP. 

Here we focus our attention on the planetary system 
WASP-44, discovered bv lAnderson et al.l lj2012l ). The system 
consists of a 0.89 Mjup planet on a 2.42-day orbit around an 
inactive G8 V star {V = 12.9, [Fe/H] = -1-0.06). In this work 
we present the first photometric follow-up since its discovery 
was announced, covering seven optical/NIR passbands. We 
refine the physical properties of the system and attempt to 
probe for radius variations in these passbands. 

Our paper is structured as follows. In §[2] we describe 
the instruments used for the observations and give some 
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details concerning the data reduction. In § [3] we illustrate 
the analysis of the data which led to the refinement of the 
orbital period and the physical properties of the WASP-44 
system. Then, we examine the variation of the planetary 
radius with wavelength. Finally, in § |4] we summarise our 
results. 



2 OBSERVATIONS AND DATA REDUCTION 

Two transits of WASP-44 b were recorded on 2011 October 
2 and 6, using the Gamma Ray Burst Optical and Near- 
Infrared Detector (GROND) instrument mounted on the 
MPgQ/ESO 2.2 m telescope at ESQ La SiUa, Chile. GROND 
is an imaging system capable of simultaneous photometric 
observations in four optical (identical to Sloan q\ r' , i' , z' ) 
and three NIR (J, H, K) passbands (iGreiner et al.ll2008l ). 
Each of the four optical channels is equipped with a back- 
illuminated 2048 X 2048 E2V CCD, with a field of view 
(FOV) of 5.4' X 5.4' at a scale of 0.158" /pixel. The three 
NIR channels use 1024 x 1024 Rockwell HAWAII- 1 arrays 
with a FOV of 10' x 10' at 0.6" /pixel. 

We applied a slight telescope defocus and obtained 
repeated integrations during both runs. Autoguiding was 
used to keep the stars on the same pixels. All images were 
digitised using the fast read-out mode (~10s) to improve 
the time sampling. On 2011/10/02 we monitored the flux 
of WASP-44 in nearly photometric conditions, and saw a 
monotonical seeing decline from 0.50" to 1.03" as airmass 
decreased from 1.08 to 1.05 and then rose to 1.54. Expo- 
sures of 20 s were used in the optical channels and stacks 
of four 4 s images were obtained in the NIR channels. On 
2011/10/06 we experienced worse atmospheric conditions, 
and the airmass changed from 2.07 to 1.05 during our obser- 
vations. Images with exposure times of 25 s were obtained in 
the optical channels, and stacks of six 4 s images in the NIR 
channels. For the optical frames we selected a fast read-out 
mode, which provided a 15 s readout. Due to the necessary 
synchronization of the optical and NIR starting times for 
each exposure, the observing cadence resulted to be of « 48 
and ~ 25 s in the optical and NIR bands respectively. 

Another full transit of WASP-44 b was observed 
through a Gunn r fllter on the night of 2011/09/02 using the 
Wide Field Camera (WFC) on the 2.5 m Isaac Newton Tele- 
scope (INT) at La Palma. WFC is an optical mosaic camera 
consisting of four thinned EEV 2kx4k CCDs, with a plate 
scale of 0.33" pixel"^. The telescope was heavily defocussed 
so the point spread function (PSF) of the target and com- 
parison stars had not more than 35 000 counts per pixel. We 
used the fast readout mode to reduce readout time down to 
about 15 s. An exposure time of 2 min per frame was used, 
and the resulting observing cadence was roughly 135 s. The 
night was photometric. We were not able to autoguide the 
telescope as the autoguider is incorporated into the WFC 
and so was also defocussed. A summary of the observational 
data is reported in Table[T] 

Reduction of the optical frames was undertaken using 
standard methods. We created master bias and flat-field 
images by median-combining sets of bias images and sky 

^ Max Planck Gesellschaft. 



Table 2. Excerpts of the light curves of WASP-44. The full 
dataset will be made available at the CDS. 



Telescope 


Filter 




BJD (TDB) 


Diff. mag. 




Error 


INT 






Gunn 


r 


2455807.576418 


0.00036 


0, 


.00041 


INT 






Gunn 


r 


2455807.579370 


0.00009 


0, 


.00042 


ESO 


2, 


.2m 


Sloan 


g' 


2455836.647443 


0.00126 


0, 


.00100 


ESO 


2, 


.2m 


Sloan 


g' 


2455836.649992 


0.00210 


0, 


.00123 


ESO 


2, 


.2m 


Sloan 


t' 


2455836.647443 


-0.00097 


0, 


.00088 


ESO 


2, 


.2m 


Sloan 


r' 


2455836.649992 


0.00079 


0, 


.00109 


ESO 


2, 


.2m 


Sloan 


i' 


2455836.647443 


0.00096 


0, 


.00109 


ESO 


2, 


.2m 


Sloan 


i' 


2455836.649992 


0.00138 


0, 


,00134 


ESO 


2, 


.2m 


Sloan 


z' 


2455836.647443 


-0.00143 


0, 


,00142 


ESO 


2, 


.2m 


Sloan 


z' 


2455836.649992 


0.00100 


0, 


,00179 



fiats, and used them to correct the science images. Aperture 
photometry was performe d using the lDl[f|/ASTROLl jf| im- 
plementation of daophot (|Stetsonlll987l : ISouthworth et al.1 
I2OO9I ). The apertures were placed manually and shifted to 
account for pointing variations, which were measured by 
cross-correlating each image against a reference image. We 
experimented with wide range of aperture sizes and retained 
those which gave photometry with the lowest scatter com- 
pared to a fitted model. The times of observation were con- 
ve rted from UTC to BJ D(TDB) using the IDL procedures 
of'East man et al.l (|2010l '). 

Differential photometry was obtained in each filter using 
between two and four comparison stars, two of which were 
enough bright to produce count rates comparable to those 
of WASP-44. All good comparison stars were combined into 
one ensemble by weighted flux summation. Slow variations 
in the apparent brightness of the reference stars, primarily 
attributable to atmospheric effects, were treated by fltting 
a polynomial to regions outside transit, whilst simultane- 
ously optimising the weights of the comparison stars. Given 
the limited number of comparison stars in the small fleld 
of view of GROND, and the shape of the slow brightness 
variations, we used a second-order polynomial versus time. 
Uncertainties introduced by this procedure were considered 
in the modelling of the light curves, as described in the next 
section. 

The NIR frames were also calibrated in a standard way, 
including dark subtraction, flat correction and sky subtrac- 
tion. The master sky images used in sky subtraction were 
created from two sets of 20-position dithering sky measure- 
ments, one before the science observation and one after. We 
performed aperture photometry on the calibrated NIR im- 
ages as well. The aperture locations were determined us- 
ing idl/find. Various combinations of aperture and annulus 
sizes were check to flnd the best photometry. We also care- 
fully made ensembles of comparison stars, and chose the 
group which showed the least deviation from the target. Af- 
ter normalising the target light curve with the composite 
reference light curve, we decorrelated the data with posi- 
tion, seeing and airmass in order to remove the correlated 

^ The acronym IDL stands for Interactive Data Language and 
is a trademark of ITT Visual Information Solutions. For further 
details see http://www.ittvis.com/ProductScrvices/IDL.aspx 
^ http://idlastro.gsfc.nasa.gov/ 
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Table 1. Log of the observations presented in this work. Af^bs is the number of observations and "Moon ilium." is the fractional illumination 
of the Moon at the midpoint of the transit. The aperture sizes are the radii of the software apertures for the star, inner sky and outer sky, 
respectively. /3 is the factor used to inflate the errorbars (see §|3ll. Transit #1 was observed using the 2.5 m INT in La Palma, whereas 
transits #2 and #3 using the MPG/ESO 2.2 m telescope in La Silla. 



Transit 


Date 




Start /end Time 
(UT) 




Exposure 
time (s) 


Filter 


Airmass 


Moon 
ilium. 


Aperture 
sizes (px) 


Scatter 
(mmag) 




/3 


1 


2011 


09 


03 


00 


28 


- 00:57 


156 


120 


Gunn r 


1.60 1.32 2.18 


31% 


25, 35, 50 


0.37 


1 


15 


2 


2011 


10 


02 


03 


21 


- 07:35 


317 


20 


Sloan g' 


1.08 1.05 -i- 1.54 


40% 


16, 35, 55 


1.23 


1 


07 


2 


2011 


10 


02 


03 


21 


- 07:35 


317 


20 


Sloan r' 


1.08 1.05 -i- 1.54 


40% 


17, 40, 60 


1.08 


1 


00 


2 


2011 


10 


02 


03 


21 


- 07:35 


317 


20 


Sloan i' 


1.08 1.05 -)• 1.54 


40% 


20, 30, 45 


1.34 


1 


00 


2 


2011 


10 


02 


03 


21 


- 07:35 


317 


20 


Sloan z' 


1.08 1.05 1.54 


40% 


12, 25, 40 


1.79 


1 


10 


2 


2011 


10 


02 


03 


21 


- 07:35 


631 


3 


J 


1.08 1.05 -i- 1.54 


40% 


6.0, 15.0, 23.0 


7.66 


1 


09 


2 


2011 


10 


02 


03 


21 


- 07:35 


631 


3 


H 


1.08 1.05 -i- 1.54 


40% 


2.5, 11.5, 21.0 


8.11 


1 


08 


2 


2011 


10 


02 


03 


21 


- 07:35 


631 


3 


K 


1.08 1.05 1.54 


40% 


3.0, 11.0, 21.0 


5.98 


1 


15 


3 


2011 


10 


06 


03 


21 


- 07:35 


317 


25 


Sloan g' 


2.07 1.05 


80% 


19, 30, 50 


1.89 


1 


06 


3 


2011 


10 


06 


03 


21 


- 07:35 


313 


25 


Sloan r' 


2.07 -)• 1.05 


80% 


22, 40, 60 


1.70 


1 


32 


3 


2011 


10 


06 


03 


21 


- 07:35 


313 


25 


Sloan i' 


2.07 -)• 1.05 


80% 


20, 40, 60 


1.79 


1 


35 


3 


2011 


10 


06 


03 


21 


- 07:35 


313 


25 


Sloan z' 


2.07 1.05 


80% 


16, 40, 60 


2.52 


1 


26 


3 


2011 


10 


06 


03 


21 


- 07:35 


623 


3 


J 


2.07 -i- 1.05 


80% 


9.0, 17.0, 25.0 


5.23 


1 


38 


3 


2011 


10 


06 


03 


21 


- 07:35 


623 


3 


H 


2.07 -i- 1.05 


80% 


6.5, 13.5, 22.0 


6.96 


1 


34 


3 


2011 


10 


06 


03 


21 


- 07:35 


623 


3 


K 


2.07 -> 1.05 


80% 


3.5, 4.5, 18.0 


7.65 


1 


04 



red noise. We extracted the optimal NIR light curves for the 
two nights according to the rras of O-C residuals and con- 
sistency of the transit depth between two nights. Details are 
given in Appendix |X] The resulting photometry is given in 
Table 2. 

Although we have performed all the possible calibra- 
tions and corrections, the scatter of the NIR data is much 
larger than that of the optical data, and the NIR light curves 
are heavily dominated by red noise. This problem is related 
to the adopted observing strategy that unfortunately did 
not allow a good SNR to be obtainecQ. 

Our analys i s also included the dataset presented by 
lAnderson erahl (|2012l ). which was obtained through a Gunn 
r filter and using EulerCam mounted on the 1.2 m Euler- 
Swiss telescope at ESO La Silla. 



3 ANALYSIS 

We have measured the physical properties of the WASP-44 
planetary system following the methodology of the Homoge- 
neous Studies project (Soutliworth 2008, 2009, 2010, 201L, 
We refer the reader to those works for a detailed de- 
scription of the approach. 

Since the absolute values of the observational errors 
from our pipeline (which come ultimately from the APER 
subroutine) were found to be underestimated, we adopt the 
standard practice of rescaling them for each dataset to give 
a reduced x oi Xv = ^- Then, in order to account for time- 
correlated errors (i.e. red noise, which can significantly affect 
ground-based data) and derive more realistic uncertainties, 
we inflated the errorbars further by multiplying the data 



* After these first observations, we changed our observing strat- 
egy adopting the defocussing technique for GROND. Thanks to 
this approach, we can now obtain scatter smaller than 2 mmag 
per observation in the NIR bands and smaller than 1 mmag in 
the optical ones, without compromising the sampling (Nikolov et 
al., in prep.). 



weigh ts by a factor /? ^ 1. The /3 approach (e.g. iPont et al.l 
,2009 ; IWinn et~aL 20071) is a wi d ely used way to assess 



red noise (e.g. iGillon et all 12 006': 'Wi nn et all |2008|. l2009l: 
Gibson et al] l2008l: iNikolov et al.. .20121 : ISouthworth et all 



2012al M. The factor /3 is a measurement of how close the 



data noise is to the Poisson approximation, and is found 
by binning the data and evaluating the ratio between the 
size of the residuals versus what would be expected if the 
data followed Poisson statistics. We evaluated values of the 
P factor for each individual transit and for groups of 10 dat- 
apoints; they are reported in Table[T] For a more exhaustive 
discussion about rescaling errorbars, see lAndrae, (,2010[ ). 



3.1 Period determination 

As a first step, we fitted each light curve individually us- 
ing the JKTEBOP code (see §[331 order to find the tran- 
sit times. Their uncertainties were estimated from Monte 
C arlo simulations. To the se we added two timings obtained 
bv lAnderson et al] (|2012l ) and four obtained by amateur as- 
tronomers, which are available on the TRESCAQ 

website 

(see Table[3]). We excluded from this analysis our NIR data 
since they are much noisier than the optical ones. All timings 
were placed on BJD(TDB) time system. The resulting mea- 
surements of transit midpoints were fitted with a straight 
line to obtain a new orbital ephemeris: 

To = BJD(TDB)2 455 434.37642(37) +2.4238133(23) x E 

(1) 

where E is the number of orbital cycles after the refer- 
ence epoch (the mid point of the first transit observed by 
lAnderson et aLll2012l ) and quantities in brackets denote the 
uncertainty in the final digit of the preceding number. The 
fit has xt ~ 2.42, so the uncertainties were increased to 
account for this. Despite this large xS, we do not note any 



5 The TRansiting Ex oplanetS and CAndidat es (TRESCA) web- 
site can be found at http://var2.astro.cz/ENj /tresca/index.php 
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Table 3. Central transit times of WASP-44 and their residuals 
versus the ephemeris derived in this work. TRESCA refers to the 
"TRansiting ExoplanetS and CAndidates" website. 



Central transit time Cycle 
BJD(TDB) - 2400000 no. 



Residual Reference 
(JD) 



55434, 
55453, 
55807, 
55814, 
55829, 
55829, 
55836, 
55836, 
55836, 
55836, 
55841, 
55841, 
55841, 
55841, 
56127, 



37637 
76639 
64374 
91655 
45489 
46151 
72905 
72979 
72900 
72928 
57719 
57757 
57684 
57769 
58624 



± 0.00040 
± 0.00042 
±0.00013 
±0.00150 
± 0.00245 
±0.00163 
± 0.00020 
± 0.00030 
± 0.00020 
±0.00015 
± 0.00035 
± 0.00046 
± 0.00028 
±0.00031 
± 0.00048 
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154 
157 
163 
163 
166 
166 
166 
166 
168 
168 
168 
168 
286 



-0.00005 
-0.00054 

0.00007 

0.00144 
-0.00310 

0.00315 
-0.00038 

0.00036 
-0.00043 
-0.00015 

0.00014 

0.00052 
-0.00021 

0.00064 
-0.00078 



Anderson et al. (_2012) 
Anderson et al. ('2012 ) 



This work (INT r) 
Evans P. (TRESCA) 
Lomoz F. (TRESCA) 
Lomoz F. (TRESCA) 
This work (GROND g') 
This work (GROND r') 
This work (GROND i') 
This work (GROND z') 
This work (GROND g') 
This work (GROND r') 
This work (GROND i') 
This work (GROND z') 
Sauer T. (TRESCA) 



systematic deviation from the predicted transit times. A plot 
of the residuals around the fit is shown in Fig.[T] Since the 
number of known observed transits of this planet is still very 
low, no conclusions can be drawn regarding the existence of 
any third body in the system. 
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Figure 4. As for F ig.[2] but for the r-filter INT and Euler 
l|Anderson et al.ll20l3 l light curves of WASP-44. 



3.2 Light curve modelling 

The light curves were modelled using the ,TKTEBOllf| code. 
The primary fitted parameters were the sum and ratio of the 
fractional radii of the star and planet, rA+^b and k = rb/rA, 
and the orbital inclination, i. The fractional radii of the com- 
ponents are defined as ta ~ Ra/cl and rb = Rh/a, where a 
is the orbital semimajor axis, and Ra and i?b are the true 
radii of the two objects. Additional parameters of the fit in- 
cluded the light level outside transit and the midpoint of the 
transit. Limb darkening (LD) was imposed using a quadratic 
law and the corresponding coefficients fixed a t values theo - 
retically predicted using model atmospheres l|Claretll200j ). 
ISouthworthI l|2008l ) performed three different treatments of 
the LD coefficients: (i) both coefficients fixed to theoreti- 
cal values; {ii) both coefficients fitted; {iii) the linear LD 
coefficient fitted and the nonlinear one fixed to its theoret- 
ically predicted value. While the second treatment is possi- 
ble only when the data are of high quality, the use of the 
third alternative has a negligible impact on the final results. 
Uncertainties were calculated using both Mont e Carlo sim- 
ulations and a residual-permutation algorithm l|Southworthl 
L2OO8), and the larger of the two values was retained for each 
output quantity. The orbital eccentricity was fixed to zero 
l|Anderson et al.ll2012h . All the datasets were solved individ- 
ually. As in the previous section, due to their large scatter 
(see Table[T]|, we did not use the NIR data for estimating 
the physical properties of the planetary system. 

Representative photometric parameters were obtained 
for each dataset and are given in Table|4] The corresponding 

^ JKTEBOP is written in FORTRAN77 and the source code is 
available at ht tp://www.astro.keele.ac.uk/~l jkt/ 



best fits are shown in Fig.[2]and|3]for the GROND data and 
in Fig.|3] for the INT and Euler data. The final photomet- 
ric parameters are the weighted mean of the values for each 
dataset. The agreement between light curves is excellent: the 
value of the agreement of the individual parameters with 
respect to the weighted mean is smaller than 0.75 for all 
except fc, where we found a modestly larger xi 1.64. We 
have found this situation to be common during our work on 
the Homogeneous Studies pa pers. Table[4]also shows a com- 
parison with the results from lAnderson etai] (I2OI2I ). which 
are in good agreement with ours but have much larger er- 
rorbars. 

3.3 Physical properties of the WASP-44 system 

Following the approach described in lSouthworthI (|2009l ). the 
estimation of the physical properties of the WASP-44 sys- 
tem w as performed making use of standard formula (e.g. 
iHildit ch 2001). We used the photometric parameters mea- 
sured in ij |3.2l and the velocity amplitude, effective temper- 
ature, and metallicity of the star {Ka = 138.8 ± 9.0kms~^, 
Teff = 5410 ± 150 K [Fe/H]= +0.06 ± 10) measured by 
lAnderson etal] (|2012| ). We interpolated within tabulations 
from theoretical stellar models to find the best agreement be- 
tween the observed and model-predicted T^s, and the mea- 
sured rA and calculated Ra/cl. This yielded the best-fitting 
mass, radius, surface gravity and mean density of the star 
(A/a, Ra, log 5a and pa) and of the planet (A/b, Rh, Qb 
and pb). We also calculated the orbital semi- major axis (a), 
planetary equilibrium temperature (Teq), Safranov number 
(O) , and the evolutionary age of the star. 

The uncertainties in the input parameters were prop- 
agated into the output physical properties using a pertur- 
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Figure 1. O-C diagram of the mid-tra nsit times of WASP-44b versus a linear ephemeris. The timings in black are from this work, and 
in grey are from I Anderson et al. 1 1|2012| '). The timings obtained by amateur astronomers are plotted using open circles. The uncertainties 
of our points have been rescaled (see text). 
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Figure 2. First set (2011/10/02) of GROND light curves of WASP-44 compared to the best JKTEBOP fits using the quadratic LD law. 
The residuals of the fits are plotted at the base of the figure, offset from zero. 



bation analysis l|Southworth et aLll2005l ). Systematic errors 
were assessed by comparin g results from five different sets 
of theoretical models (see ISouthworthI |2010| ). The sets of 
physical properties found using each set of stellar models 
is shown in Tableland the final physical properties of the 
WASP-44 system are given in Table[6l The results obtained 
by [Anderson et al., (,2012|) are less precise but are in good 



agreement with our own. In particular, we improved the 
measurement precisions of the radii of both the planet and 
the parent star, by factors of 4 and 3 respectively. The radius 
of WASP-44 b that we found (1.002 ± 0.033 ± 0.018 Rjup) is 
smaller than that measured in the discovery paper and does 
not match the pre dicted radii of corcless hot- Jupiter plan- 
ets as estimated bv lFortnev et alj (|2007i ) . According to their 
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Figure 3. Second set (2011/10/06) of GROND light curves of WASP-44 compared to tlie best jktebop fits using ttie quadratic LD law. 
The residuals of the fits are plotted at the base of the figure, offset from zero. 



Table 4. Parameters of the fits to the ten light curves of WASP-44. The final parameters are the vireighted mean of the result for the 
GROND, INT and Euler light curves. 



Source 




k 


i° 


rA 


rh 


GROND g'-band #1 
GROND g'-band #2 
GROND r'-band #1 
GROND r'-band #2 
GROND i'-band #1 
GROND i'-band #2 
GROND z'-band #1 
GROND z'-band #2 
INT r-band 
Euler r-band 


0.1368 ±0.0056 

0.139 ±0.010 
0.1268 ±0.0045 
0.1284 ±0.0080 
0.1363 ±0.0051 
0.1303 ±0.0087 
0.1306 ±0.0075 

0.123 ±0.010 
0.1285 ±0.0029 
0.141 ±0.011 


0.1228 ±0.0023 
0.1188 ±0.0028 
0.1175 ±0.0014 
0.1190 ±0.0029 
0.1205 ±0.0016 
0.1204 ±0.0029 
0.1153 ±0.0023 
0.1134 ±0.0034 
0.1196 ±0.0013 
0.1197 ±0.0039 


86.2 ±0.5 

86.2 ± 1.0 

87.3 ±0.6 
86.9 ±0.8 

86.3 ±0.5 

86.7 ±0.9 

86.8 ±0.8 

87.4 ± 1.5 
86.6 ±0.3 

85.9 ± 1.0 


0.1218 ± 0.0020 
0.1240 ±0.0090 
0.1135 ± 0.0039 
0.1147 ±0.0068 
0.1216 ± 0.0044 
0.1163 ± 0.0076 
0.1171 ± 0.0065 
0.1102 ± 0.0089 
0.1148 ± 0.0026 
0.1261 ± 0.0100 


0.01496 ± 0.00086 
0.0147 ±0.0014 

0.01333 ± 0.00060 
0.0136 ± 0.0011 
0.0147 ± 0.00068 
0.0140 ±0.0011 

0.01349 ± 0.00095 
0.0125 ± 0.0013 

0.01373 ± 0.00042 
0.0151 ± 0.0016 


Final results 






86.59 ±0.18 


0.1168 ±0.0016 


0.0139 ±0.0002 


Anderson et al. (2012) 


0.1398 


0.1260 ±0.0030 




0.1242 ± 0.0102 


0.0157 
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Table 5. Derived physical properties of WASP-44 from using each of five different theoretical stellar models. In each case gi, = 21.5 ± 
1.6ms-2, PA = 1.414 ± 0.058 P0 and Toq = 1304 ±37K. 





This work 


This work 


This work 


This work 


This work 




(Claret models) 


(Y2 models) 


(Teramo models) 


(VRSS models) 


(DSEP models) 


i^Tb (kms-i) 


156.1 ±4.1 


151.7 ± 1.0 


152.6 ±4.2 


152.7 ±4.1 


153.7 ±3.2 


Ma (Mq) 


0.968 ± 0.077 


0.888 ±0.018 


0.904 ± 0.074 


0.907 ±0.073 


0.924 ±0.058 


Ra {%) 


0.881 ±0.025 


0.856 ± 0.016 


0.862 ± 0.024 


0.862 ±0.024 


0.865 ±0.020 


logSA (cgs) 


4.534 ±0.018 


4.521 ± 0.010 


4.524 ±0.018 


4.524 ±0.018 


4.530 ±0.016 


Mb (Mjup) 


0.901 ±0.075 


0.851 ± 0.056 


0.861 ± 0.073 


0.863 ± 0.073 


0.874 ± 0.067 


Rh (^jup) 


1.020 ±0.033 


0.992 ± 0.019 


0.997 ±0.033 


0.998 ± 0.032 


0.994 ±0.027 


Ph (Pjup) 


0.793 ±0.070 


0.816 ± 0.069 


0.812 ± 0.072 


0.811 ±0.072 


0.831 ±0.071 


e 


0.0640 ± 0.0046 


0.0658 ± 0.0045 


0.0654 ± 0.0048 


0.0654 ± 0.0047 


0.0656 ± 0.0046 


a (AU) 


0.03508 ± 0.00093 


0.03409 ± 0.00023 


0.03429 ± 0.00093 


0.03432 ± 0.00092 


0.03454 ± 0.00072 



Table 6. Final physical properties of the WASP-44 system. The 
first errorbar for each parameter is the statistical error, which 
stems from the measured spectroscopic and photometric parame- 
ters. The second errorbar is the systematic error arising from the 
use of theoretical stellar models, and is given only for those pa- 
rame ters which have a depe ndence on stellar theory. The results 
from [Anderson et al.l l|2012l 'l are included for comparison. 



This work (final) 



Anderson et al. 



Ma (Mq) 0.917 ± 0.077 ± 0.051 

Ra (Rq) 0.865 ±0.025 ±0.016 

log gA (cgs) 4.526 ± 0.018 ± 0.008 

PA (pq) 1.414 ±0.058 

(Mjup) 0.869 ± 0.075 ± 0.032 

Rh (Rjup) 1.002 ±0.033 ±0.018 

£/b (ms-2) 21.5 ± 1.6 

Pb (Pjup) 0.808 ± 0.072 ± 0.015 

Teq (K) 1304 ± 37 

e 0.0652 ± 0.0048 ± 0.0012 

a (AU) 0.03445 ± 0.00093 ± 0.00063 

Age (Gyr) ^-itl^^tH 



0.951 ±0.034 

u.azi _o.074 
4 401 +0.068 

1 in+0.32 

0.889 ± 0.062 
1.14 ±0.11 

7+3.4 
-3.0 
,+0.23 
-0.15 
1343 ± 64 



15.7: 

0.61^ 



0.03473 ± 0.00041 



tablefQ, a 0.875 Mj planet with 50 A/Earth core at 0.045 AU 
from the Sun (age 3.16 Gyr), is 0.991 -R,7up. On the contrary, 
the expected radius for a core-free planet with the same 
properties is 1.119 J?,iup. 

3.4 Variation of planetary radius with wavelength 

As an additional possibility for the GROND data, we made 
an attempt to investigate the possible variation of the radius 
of WASP-44 b with wavelength. This is quite a difficult task 
because the relative faintness of the host star {V ~ 12.9) 
makes this system less than ideal for optical photometric 
studies, especially from the ground. Moreover, as already 
mentioned, our m- focus observing strategy has proved to 
be ill-suited to obtaining high-precision photometry, partic- 
ularly in the NIR bands. Finally, by using broad-band filters 
we are forced to average our measurements of the planet ra- 
dius in each band over a quite larger range of the wavelength 
(the best case in the i' filter where its width is "only" 100 

^ The tables, interpolate d from the giant plan et thermal evolu- 
tion models described in iFortnev et al.l | |2007|) , are available at 
|http: / /www.ucolick.org/ ~ jfortney/models.htm 



nm). For these reasons, the results of this section should be 
used with care. 

The GROND instrument was conceived only for the 
follow-up of gamma-ray bursts, and was not designed to al- 
low filter changes for individual observing sequences. We 
were therefore unable to use filters with narrow pass- 
bands covering specific wavelength ranges. We note that 
the four optica l bands of GROND were already used by 
Ide Mooii etHI l|2012l ) to investigate the atmosphere of 
GJ 1214 b, a 6.55 M® transiting planet. Here we try for the 
first time to use all the seven bands, exploiting the full po- 
tential of GROND. 

Fo llowi ng the same strategy used bv ISouthworth et al.l 
l|2012bl ) and lMancini et al] l|2013l ). we proceeded as follows. 
The two GROND datasets were combined by phase and ac- 
cording to passband, and then fitted with all parameters 
fixed to the final values given in Table]!] with the exception 
of k. The LD coefficients were fixed to theoretical values. 
The er rors were estimat ed by a residual-permutation algo- 
rithm (|Southworthll2008l ). This approach removes sources of 
uncertainty common to all datasets, allowing us to maximise 
the accuracy of estimations of the fractional planetary ra- 
dius rb = Rh /cL as a function of wavelength and with relative 
errorbars only. The results are displayed in Fig.O where the 
points show the data, the vertical bars represent the rela- 
tive errors in the measurements and the horizontal bars in- 
dicate the full widths at half maximum transmission of the 
passbands used. The trasmission curves are also reported 
for completeness. As expected, the uncertainties in the NIR 
bands are larger, due to the larger scatter and systematic 
features in the hght curves. 

Inspection of Fig. [5] shows that our measurements are 
unfortunately not sufficiently accurate to claim a clear vari- 
ation of rb along the seven passbands. Using the plane- 
tary system values reported in Table[6l we computed a one- 
dimensional model atmospher e of WASP-44 b, us i ng th e 
atmosphere code described in iFortnev et all (|2005l . I2OO8I ). 
The fully non- gray model uses the che mical equilibrium 
abundances of iLodders &: FeglevI ||2002|) and the opacity 

J (12 



database described in Freedman et al. riiooi). The atmo 



spheric pressure-temperature profile simulates planet-wide 
average conditions. We furthermore computed the transmis- 
sion spectrum of the model using the methods described in 
IFortnev et al] (|2010l ). 

At optical wavelengths, rb appears to be slightly larger 
in r' and i' than g' and z' . This is as expected for a planet 
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Figure 5. Variation of the fractional planetary radius rb = Rh/o, with wavelength. The points shown in the plot are from the MPG/ESO 
2.2 m telescope. The vertical bars represent the errors in the measurements and the horizontal bars show the full widths at half maximum 
transmission of the passbands used. The solid, blue continuous line is the calculated synthetic spectrum based on a theoretical model of 
the atmosphere of WASP-44 b. Red boxes indicate the predicted values for this model integrated over the passbands of the observations. 
Transmission curves of the GROND filters are shown at the bottom of the figure. 



with an equilibrium temperature Toq ~ 1300 K, due to the 
larger predicted radius around the Na and K lines. Although 
a similar effect is seen in the model (Figure[6ll , our data are 
not sufficiently precise to reject the default hypothesis that 
Tb is constant with wavelength. The radius of WASP-44 b 
is in agreement with the theoretical model in the J band, 
but is larger in the H and K bands. We found that the 
in r' is 9.5% and 4.7% smaller than that for the H and K 
band respectively, which correspond to a difference in units 
of atmospheric pressure scale height (H) of roughly 23H 
and lOH, respectively. The first difference is very large, and 
would involve an extreme opacity in H . 



A more plausible alternative is that the NIR photome- 
try is affected by systematics or correlated noise, which was 
not possible to eliminate. This is also suggested by the poor 
SNR observed in the data, which is reflected in the higher 
scatter of the points in the J, H, K light curves and the dif- 
ferent shapes of the fit in the same band between the first 
and the second observed transits (Figs.[5]and|3}. We deduce 
that systematic errors in the NIR dominate the error budget 
at these wavelengths. 



4 SUMMARY AND DISCUSSION 

In this work we have presented the first photometric follow- 
up of the transiting extrasolar planetary system WASP-44. 
One transit was observed using WFC on the INT, and two 
transits were observed in seven passbands simultaneously 
using GROND on the MPG/ESO 2.2 m telescope. Using 
these new datasets plus ano ther one taken from the liter- 
ature ([Anderson et al.l |2012| ) , we have obtained improved 
measurements of the orbital ephemerides and physical prop- 
erties of the system. In particular, we found that the radius 
of WASP-44 b is smaller than previously found, and is differ- 
ent (at 2a confidence le vel) from that theoretically expected 
for a free-core planet (|Fortnev et al ] 12007. ) . This suggests 
that WASP-44 b is a "heavy element rich" planet, with the 
heavy elements being in a distinct core or mixed within the 
H/He envelope. Such a result is important in order to draw 
an accurate mass-radius plot for TEPs, which provides key 
diagnostic for theoretical works that look to infer the bulk 
composition of the giant planets and distinguish them (in 
some cases) from brown dwarfs. 

Although GROND was not designed for follow-up stud- 
ies of planetary transits, it is one of the very few imaging sys- 
tems worldwide able to perform simultaneous photometric 
multi-band observations. The later fact makes GROND an 
efficient and very useful instrument to detect transit anoma- 
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Figure 6. As in Fig. [5] but zoomed in on the optical wavelengths. In this range, the synthetic transmission spectrum is dominated by 
gaseous Na and K. 



lies or to probe the atmosphere of TEPs. Indeed, atomic 
and molecular absorption as well as scattering processes 
may result in detectable rad i us va r iations as a functio n 
of wavelength ("Portnev et ajj l2007l : iBatygin etHI |2009| ). 
Ide Mooii eta l. (2012) already used GROND in the case of 
GJ 1214 b, but they restricted their analysis only to the opti- 
cal bands. Here, for the first time, we used all the seven band 
of this powerful instrument, to search for variations of the ra- 
dius of WASP-44 b with wavelength. We measured the plan- 
etary radius in the seven GROND passbands, correspond- 
ing to a wavelength coverage of 370-2440 nm, and compared 
them to a synthetic spectrum based on an isothermal model 
atmosphere in chemical equilibrium. The model reproduces 
the radius values in the optical, in agreement with the equi- 
librium temperature of WASP-44 b but, our data lack the 
precision to rule out the possibility that the planetary ra- 
dius does not vary with wavelength. It does not predict the 
radius values found in the H and K bands, which differ 
by roughly 23 and 10 atmospheric pressure scale heights to 
those in the other bands. This phenomenon is most likely 
due to the comparatively low quality of the NIR data. We 
have worked to reduce the sources of noise in the subsequent 
uses of the GROND instrument. More enticing results were 
obtained for other TEPs and will be shown in forthcoming 
papers. 
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APPENDIX A: GROND-NIR DATA 
REDUCTION 

The calibration of the GROND NIR data is a little more 
complicated than the optical data, because of the presence 
of electronic odd-even read-out pattern along the y-axis. In 
order to remove this pattern, we smooth each image and 
compare it with the unsmoothed ons after the master dark 
has been subtracted. The amplitudes of the read-out pattern 
are then determined by comparing the median level of each 
column to the overall median level. Each column is shifted 
back to the overall median level. The master skyflats are 
divided out from the science images after the removal of 
read-out pattern. 

Since our observations are in staring mode, sky back- 
ground is an important contribution in the images. With the 
20-position dithering sky measurements right before and af- 
ter the science time series, we could construct a sky emis- 
sion model for each science image. These sky images are 
calibrated in the same way as the science images. We mask 
out all the stars in the calibrated sky images, and normalise 
each image with their own overall median levels. All the 
normalised skies are then median combined, and normalised 
again. For each science image, we scale the pre- and post- 
science sky models to their background individually. The 
final sky model is constructed by weighted combining them 
according to the inverse square of fitted . This sky model 
is then subtracted from corresponding science image. After 



sky subtraction, we still see remnant structure, which is ex- 
pected due to the variation of the sky during the observing 
sequence. 

After these calibrations have been applied, wc perform 
aperture photometry on WASP-44 as well as three nearby 
comparison stars. The locations of each star are determined 
by idl/find, which calculates the centroids by fitting Gaus- 
sians to the marginal x and y distributions. In order to find 
the optimal photometry, we lay 30 apertures on all the stars 
in step of 0.5 pixel, each with ten annuli in step of 1 pixel. 
This produces 300 datasets with different aperture settings. 
In this way, the above mentioned sky subtraction remnant 
effect is accounted for in the skj' annulus adopted in the 
aperture photometry. We divide the flux of each star with its 
own median value. A composite reference light curve is con- 
structed by median-combining fluxes of ensembles of com- 
parison stars, which are required to show the least deviation 
from the target. We normalixe the target light curve with 
this composite reference light curve. 

Since the normalised light curve still shows strong red 
noise correlated with positions, seeing and airmass, instead 
of selecting the optimal photometry directly from rms of 
out-of-transit baseline, wc fit the whole light curve simul- 
taneously with the theoretical light curve model multi- 
plied with a baseline correction function. This correction 
function consists of positions/time in quadratic form, and 
FWHMs/airmass in linear form. We tried to select the op- 
timal light curves by choosing the ones with the least rms 
of O-C flux residuals among all the datasets. However, the 
transit depths of the optimal light curves of the same band 
from different nights are very deviant from each other. We 
note that light curves of the same band obtained using dif- 
ferent aperture and annulus sizes could exhibit similar least 
rms, however produce different transit depths. Therefore, 
for each NIR band, we choose the two light curves with the 
least rms and with consistent transit depth between the two 
nights as our final optimal photometric results. 
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